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Abstract 

The orthopositronium decay to three photons is studied in the phase-space region where one 
of the photons has an energy comparable to the relative three- momentum of the e~^e~ system 
(u ~ ma). The NRQED computation in this regime shows that the dominant contribution arises 
from distances ~ l/y^mtJ, which allows to treat the Coulomb interaction perturbatively. The small- 
photon energy expansion of the 1-loop decay spectrum from full QED yields the same result as the 
effective theory. By doing the threshold expansion of the 1-loop QED amplitude we confirm that 
the leading term arises from a loop-momentum region where ~ q^/m ~ to. This corresponds to a 
new non-relativistic loop-momentum region, which has to be taken into account for the description 
of a non-relativistic particle-antiparticle system that decays through soft photon emission. 
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I. INTRODUCTION 



Effective field theories (EFT's) and asymptotic expansions have become standard tools 
in the description of the radiative decay of bound states. Their success hes on their ability 
to select the relevant dynamics in the different kinematic regions defined by the relative size 
of the radiated photon energy as compared to the bound state mass. Methods based on 
the use of operator product expansions and effective Lagrangians do not only provide more 
efficient ways to perform the computations but also allow to extract physical interpretations 
which may go unnoticed in a full theory approach. 

Apart from being of physical relevance in itself, the Positronium system provides a test- 
ing ground for the EFT concepts and techniques devised in the study of radiative decay 
amplitudes, that could eventually be applied to the description of the more intricate radia- 
tive decays of quarkonia. In this paper we illustrate which are the characteristic features of 
the radiative decay spectrum when the energy of the photon lies in the soft-energy region 
{uj ~ ma) through the study of the 87 annihilation amplitude of the Positronium spin triplet 
ground state (orthopositronium: o-Ps). The soft-energy region is accessible in the o-Ps— > 87 
decay because one of the final state photons can have an arbitrarily small energy, the other 
two being hard photons with energy < m. The process can be then viewed as the radiative 
version of the o-Ps-^ 27 decay. 

In a previous paper [1] the o-Ps differential decay spectrum was calculated in the region 
where the photon energy is comparable to the Positronium binding energy, uu ~ ma^. Bind- 
ing effects were included for the first time in the Ps structure function using non-relativistic 
QED (NRQED) [2], and it was found that they are essential to reach agreement with the 
Low's theorem prediction for the uj ^ behaviour of this decay [1, 3]. Only the leading term 
in the multipole expansion of the radiated photon field (dipole approximation) was required 
for the calculation, in accordance with the NRQED velocity counting rules for photons with 
cu ~ ma^ [4]. Based on general considerations about the power-counting of higher-order 
multipoles in the non-relativistic description of the o-Ps decay process, it was then argued 
by Voloshin [5] that the dipole approximation used for photon energies uj ~ ma^ should 
also hold for a description of the o-Ps photon spectrum in the whole range u <ti m, thus 
enlarging the validity region of the formula given in Ref . [1] . In particular, it was shown by 
Voloshin that the expansion of the above-mentioned formula in the region ma^ <^ ou <^ m 
is actually a series in a^fmfuj ~ ^/a rather than in integer powers of a. 

The origin of such unnatural expansion has to be traced back to the fact that the 
main contribution in the EFT calculation arises from distances in the e^e~ system of 
order r ~ Xj \/mio^ that are much smaller than the size of the Ps atom r ~ 1/ma. 
Translated into momentum space, this implies that there is a kinematic region where the 
relative 3-momentum of the non-relativistic e+e^ pair obeys the non-rclativistic relation 

~ p^/m ~ u). This scaling does not correspond to any of the known modes that have 
been identified in the common apphcations of non-relativistic EFT's, and should be com- 
pared to the scaling of the familiar potential modes, ~ p^/m ~ ma^, characteristic of 
the heavy particles that form the bound state, and that of the soft modes, p° ~ u; ^ p^/m. 

In view of these novel features, it is worth questioning if the non-relativistic EFT frame- 
work devised to describe the o-Ps decay for ui ~ ma^ in Ref. [1] also provides the proper 
expansion when extended to the whole uj <^ m region, as advocated by Voloshin. It is 
the main purpose of this work to confirm that this is indeed the case by comparing the 
analytic results for the EFT and full QED spectrum in the lu ~ ma region, where, as we 
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argue, the usual perturbative QED expansion can already be applied. Such comparison 
has become feasible after a recent evaluation of the 1-loop QED o-Ps— 87 phase-space dis- 
tribution [6], that we use to determine the QED spectrum for soft photon energies. The 
agreement between both computations shows that the soft-energy region provides a regime 
where the EFT and the perturbative QED calculations can have a smooth matching. More- 
over, we shall confirm by explicitly computing the o-Ps— 87 decay spectrum in NRQED 
beyond the dipole approximation that the {uj/mY suppression advocated by Voloshin for 
the higher-order multipoles is of application in the soft-energy region. 

It is also the aim of this work to show that the momentum scale po p'^/m '-^ cu that 
rules the behaviour of the corrections to the NRQED o-Ps decay amplitude, corresponds to 
a new loop-momentum region that has to be taken into account for a successful application 
of the threshold expansion method [7] to QED and QCD loop diagrams involving a particle- 
antiparticle system decaying through the emission of a soft photon. In the conventional loop 
expansion, soft photon radiation from a heavy particle-antiparticle system introduces a soft 
energy component ui ~ ma in the zero component of the heavy particle momenta. Non- 
relativistic poles in massive propagators of the form [p^ — /2m — can be thus found in 
these loops, giving rise to a contribution from the loop-momentum region po ~ p^/m ~ cu. It 
is a well-known fact that momentum regions that have not been considered previously may 
become relevant for some kinematic configurations or at higher loop order in the asymptotic 
expansion of integrals based on the "method of regions" [22], and looking for missing regions 
is a required task in order to check the robustness of the method. Likewise, it is of conceptual 
importance to understand the role of this new momentum scale in the EFT frameworks that 
describe the radiative decays of heavy particle-antiparticle systems in the soft-energy region, 
which we discuss in this paper. 

The outline of the paper is as follows. In Sec. II we detail the EFT calculation of the o- 
Ps— > 87 decay amplitude without multipole expanding the electromagnetic interaction. We 
also show how relativistic corrections that contribute to the decay amplitude in the whole 
photon energy range uo <^m can be accommodated in the EFT framework. In particular, the 
recoil of the intermediate virtual e+e^ pair, which yields also SJiuj/m correction, is included 
in our calculation. In Sec. Ill we perform the a; ^ limit of the 1-loop QED spectrum 
directly from the 1-loop phase-space distribution computed recently by Adkins [6]. The 
1-loop QED amplitude in this regime is also obtained with the asymptotic expansion of the 
graphs near threshold in Sec. IV. A discussion on the EFT description of the new loop- 
momentum region found is postponed to the end of the latter section. Finally in Sec. V 
we suggest the use of the threshold expansion within the NRQCD factorization framework 
to compute the short-distance coefficients in quarkonium radiative decays at soft energies. 
The Appendix A collects the exact formulas for the multipole calculation of the NRQED 
spectrum, while Appendix B shows another example of the asymptotic expansion method 
applied to a radiative amplitude. 

II. NRQED COMPUTATION OF THE DECAY SPECTRUM FOR u; <C m 

The NRQED framework used in Ref. [1] provides a systematic way to compute bound 
state effects in the the o-Ps decay spectrum when the photon energy is much smaller than 
the electron mass. Contrary to the usual relativistic approach where the 8-photon annihi- 
lation is considered to take place at very short distances as compared to the range of the 
electromagnetic binding force between e'^e" , the non-relativistic description takes into ac- 
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count that there is a long-distance part in the o-Ps decay process when one the photons in 
the final state is not hard {uj ^ m). In the latter case, the decay proceeds in two steps: 
the low-energy photon is first radiated from the bound state, making a transition from the 
C-odd ground state o-Ps {^Si) to a C-even Positronium state, which subsequently decays 
into two photons (see the EFT graph of Fig. la). The decay amplitude of the intermediate 
C-even e~^e~ state has no long-distance contribution, since both photons must be hard. The 
emission of the low energy photon from o-Ps is described by the Coulomb Hamiltonian of 
the e'^e~ system in interaction with a quantized electromagnetic field: 

H — Ho + Hint : 

P2 p2 a 

Ho = -; h He , He = , 

4m m r 

£ 6 

Hint = Pi ■ A(xi) H p2 ■ A(x2) - fxcr^- B(xi) - cr^ ■ B(x2) , (1) 

m m 

where we have used the center of mass variables (r = |x|) 

P = (Pi - P2)/2 , x = xi - X2 , P = pi + p2 , X = ^^l±i^ ^ 

with the subindices 1, 2 referring to the electron and the positron, respectively, and cr^, cr^ 
being the Pauli matrices acting on the electron and positron spinors {/i — e/2m). The 
terms shown in the Hamiltonian i^int are the leading ones in the non-relativistic expansion. 
Relativistic effects can be included through higher dimensional operators suppressed by 
powers of 1/m (see e.g. [8]), but will not be needed for the purposes of this work. 

Since the computation of the photon spectrum in Ref. [1] was intended for photon energies 
comparable to the Positronium binding energy, u ~ ma^, the interaction Hamiltonian was 
used in the dipole approximation limit. This approximation amounts to evaluating the 
vector potential A in the center of mass of the Positronium system [i.e. at x = 0), which is 
fully legitimate for radiated photons with wavelengths much larger than the characteristic 
size of the Positronium atom (a = 2ao = 2 /ma). Higher order multipoles arise as a Taylor 
series in the relative coordinate x and are suppressed by powers oiuo/a under the assumption 
that the relevant amplitudes one has to compute involve integrations to spatial extents of 
order ~ a. For photons with larger energies, ui ^ ma^, this premise will invalidate the use 
of the multipole expansion. 

However, in the case of the o-Ps system that undergoes a radiative transition before 
decaying, and as it has been properly pointed out by Voloshin [5], the scale [ma)^^ does not 
constraint the maximum photon energy for which we can apply the multipole expansion. 
The reason is that after the soft photon is radiated we have to consider all possible e+e~ 
intermediate states with the right quantum numbers. The propagation of these e'^e" states 
is described by the Green's function obtained from the Hamiltonian Hq, 

Hc + -]G{^,y,K)^S{^-y) , (2) 
m J 

at energy 

-K^ /m^ Eq- u - /Am, (3) 



4 



(sec Eq. (8) below). Since the annihilation of the intermediate e'^e~ pair into two photons 
takes place at small distances, the amplitude for the full process is thus given by a convolution 
of the o-Ps ground state and the Green's function. In coordinate space, the Green's function 
has a characteristic size ruled by the exponential factor exp(— /tr). We note that for soft 
photons of energy ma^ -C a; -C m, this exponential factor constrains the relevant integration 
region to distances of order ~ y/muj, much smaller than the spatial extent of the 
initial Positronium atom. Therefore, the characteristic distance that enters in the multipole 
expansion of the soft-photon emission in the o-Ps— > 87 decay amplitude is determined by 
the falloff of the intermediate Green's function rather than by the si ze of the Ps atom, so the 
series of multipoles is actually an expansion in powers of cor ~ ^/uJ/m which can be used 
as long as a; <C m. Note also that for the intermediate pair, an iterative computation of 
the Coulomb Green's function, G = + G^VqG^ + . . . , with G^ the free Green's function, 
shows that adding a Coulomb exchange generates a term 



/ 



dryiVcG' ^aW— <1, for ma^ <^ uj <^ m , 



as G^ ~ mexp(— fi;r)/r. Therefore the Coulomb interaction can be treated perturbatively 
in the virtual e~^e~ system after soft photon radiation. 

It is illustrative to show how the latter is realized in the actual computation of the o-Ps 
decay spectrum without multipole expanding the electromagnetic potential of the radiated 
photon. For this purpose let us now consider the NRQED calculation of Ref. [1] and keep 
the full Xj dependence of the various terms in the Hamiltonian (1). We start with the 
p • A electric amplitTidc^. Using time-ordered perturbation theory (TOPT), the amplitude 
corresponding to the graph of Fig. la reads 

^ ^ (0| \n, P) ^ (n, P| {pi • A(xO - • A{^,)} |o-Ps) 

n,P 

The sum above extends over all discrete and continuum states \n, P) of the spectrum of the 
unperturbed Hamiltonian Hq, that can be written as the direct product of a plane wave with 
the cm. momentum P and a wavefunction describing the relative motion. In configuration 
space, (xiX2|n, P) = e'^'-^ Vnlx), where '?/'„(x) is an eigenstate of He with eigenvalue En, 
and we have omitted the spin wavefunction. The energy of this intermediate (virtual) e'^e~ 
state is En,p = En + P^/4m, and Eq = —ma^ /A is the o-Ps ground state energy. 

The first matrix element in the r.h.s of Eq. (4) gives the short-distance part of the 
annihilation. The quantity A^'^'^'> is the e^e" — > 27 amplitude in QED calculated as an 
expansion in the momenta of the leptons: 

A(27) = £p (Wo + Wi ■ p + W2 • P) 0p + 0(p2) , (5) 

where 0p, arc the Pauli spinor fields that annihilate the electron and positron with 
relative momentum with respect the cm. p and — p, respectively. Only the first order 
in the momentum expansion will be needed in this work. Expressions for Wq, Wi can be 



^ We refer to the amplitudes as of 'electric' type if there is a change in the parity between the atomic states 
(of 'magnetic' type otherwise). 
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(a) 



o-Ps 





0-Ps 



FIG. 1: NRQED graphs for the o-Ps^ 87 decay. The aUowed intermediate states n correspond 
to those of the spectrum of Hq. The zig-zag hne represents the low-energy photon. The black 
dot in (a) denotes either a p • A or a cr • B interaction. Graph (6) represents the radiation of the 
low-energy photon directly from the 27 annihilation vertex. 



found in Ref. [1], while W2 = (Wi + (ei • e2)cr)/2, with ei, €2 the polarizations of the 
outgoing hard photons. The W2 ■ P term in Eq. (5) is an operator which depends on the 
total momentum of the e"'"e~ pair. In principle, we need to consider the latter because in the 
o-Ps decay the intermediate e"'"e~ pair can recoil after the soft photon is radiated. Inserting 
the operators in Eq. (5) into the short-distance amplitude one gets: 



(0| |n,P> = Tr [(Wo^n(O) - z Wi ■ {V My)}y=o + W2 ■ P V^n(O)) 

The trace above is taken over spinor indices, and is the spin wavefunction of the state 
1 c 



(6) 



\n, P), which can be in a spin- state {(px^ I/V2), 

or in a spin-1 state with polarization 



vector £ s* ■ cr/V2). 

The matrix element for the radiation of the soft photon with three-momentum k and 
polarization vector e (we pick a gauge where photon polarizations are purely transverse, 
Ci ■ ki = 0, e° = 0) can be calculated easily: 

(n, P| pi • A(xi) - p2 ■ A(x2) |o-Ps) = -I (277)35(3) (p + k) (£ ■ s*) 



with the initial o-Ps state (xiX2|o-Ps) = 
state wave function in position space 

V'o(x) 



X j d^xij*ix) [e-'^-''/^ + e^^''^'^) e ■ VM^) , (7) 
i'oi^x), once we set X = 0. ipoi^) is the ground 



-x/a 



and £0 is the spin-1 polarization vector. The photon wavefunctions e^**^'^/^ arise from the 
emission of the photon from the electron or from the positron line, respectively. The p • A 
interaction can change the orbital angular momentum but not the spin of the e~^e~ pair, 
which thus remains in a spin-1 state. The delta function in Eq. (7) sets the recoil momenta 
of the intermediate e+e~ system to — k, so that En,F = En + uj^/Am is its energy ^. Taking 
into account the sum over polarizations. 



E 



6. 



^ Let us recall that in TOFT three-momentum is conserved at the vertices and virtual states are always 
on- shell. 
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we can write 



I n ^-'•^ + ^-^0 Jy=0 

+ ^ (x^ W. . k0).„ / ci^x Ie #X^^1 (^^^'■'''' + ^^'■^'') ^ • V^o(x) . (8) 

m J \^ ^ -C/n.k + o; — £/o J 

The term Wq in Eq. (6) has disappeared from the electric amphtude because it only con- 
tributes to the 27 annihilation of a spin-singlet (parapositronium: p-Ps) state. Likewise, the 
second line in Eq. (8) vanishes because the vector integration can only be proportional to the 
photon momentum k and k ■ e = 0. Therefore we are left only with the term proportional 
to Wi, where the term between brackets can be written in terms of the Coulomb Green's 
function at energy —K^/m — Eq — uj — uj"^ /Am, 



(9) 



The derivative acting on the Green's function picks out the i — 1 component of its partial 
wave decomposition in y = 0, 

00 

G (x, y,k) = Y, {21 + 1) {xy)' Pe (x • y/xy) Ge {x, y, k) . (10) 

(See Appendix C of Ref. [1] for explicit expressions of the partial waves Ge). Therefore only 
P-wave states contribute in the sum over virtual fluctuations in Eq. (4). Using Vt/'o(x) = 
— ■0o(x) xyar, we can rewrite the electric amplitude as 

M = -e^o(O)(x^W10),„e^D^^(k), (11) 

with 

i^^^(k) = - A / rf3^ ^ (g-ik.x/2 ^ gik.x/2N ^^^Q^ ^. ^) ^-r/a (^3) 

ma J r ^ ' 

The general structure of £)*-' (k) reads 

L>^^'(k) = V{u) 5'^ + V{uj)^ . (13) 

The term proportional to 'D{uo) vanishes when contracted with the photon polarization e. 
The coefficient V^u) can be projected out by contracting with the tensor 5"^^ — k*k-^/a;^, and 
after performing the angular integration, we obtain: 



P(c.) = -A /drr3G,(0,r;.)e-^/» [ A sin (^^) - J:^ cos (^) 



(14) 



Note that all recoil effects up to this order are properly accounted for by the a;^/4m term in 
the parameter k, that enters the Green's function, and are naturally suppressed with respect 
leading terms in the non-relativistic regime. 
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In addition, gauge invariance requires that the momentum operators in the short-distance 
amphtude A^'^'^^ are replaced by the covariant derivatives pj =F eA(xj) [1]. This generates 
a purely local term that has to be added to the electric amplitude calculated with the 
interaction Hamiltonian of Eq. (1). The corresponding graph is depicted in Fig. lb. This 
contribution reads 

(0| |o-Ps> = -e Vo(0) (x^ Wi • e0),„ , (15) 
and the total electric amplitude is then equal to 

Me = -eV^o(O)(x^Wi-e0),„(l + I)(a;)). (16) 

The expression for T>{u;) in Eq. (13) allows us to identify the typical distances that control 
the falloff of the integrand. The £ — 1 Green's function, 

TD K*^ TTirv 1 

Gi(0, r; k) = --e— r(2 - u) U{2 - u, 4, 2/.r) , u^ — ^-, (17) 
OTT Zk aK 

contains the exponential factor e"'"'. On the other hand, we have the factor e~^/" from 
the wavefunction of the initial state. Therefore, distances in the integrand for the electric 
amplitude are constrained by the exponential factor 

exp|-(^Av+ l/a^} , (18) 

and the qualitative behaviour of T>{lj) depends on the scaling of the soft photon energy: 

• a; ~ ma"^ {k ~ ma): both exponential factors arc of the same size, constraining the 
integration variable to r ~ a. The multipole expansion is applicable here as the product 
ujr ^ a, and it is realized in the full result (14) by expanding the trigonometric functions 
between brackets: 



24 . /a;r\ 12 fujr\ (a;r) 



2 



, sin — - — — cos — = 1 - + 0(a;V) , (19) 

and thus 

V{uj) = de{uj) + S^\uj) + . . . 

= 4(a;)(l + 0(a;Va2)) , (20) 

where the leading term is given by 

de{uj) = -^J drr^Gi{Q,r;n)e-'/\ (21) 

The quantity 

1 + de{ou) = ae{u}) 

gives the electric amplitude in the dipole hmit, which was calculated in Ref. [1] (Eq. (60) 
therein) using an interaction Hamiltonian written in a gauge invariant form in terms of the 
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electric field ^. A formula for ac{uj) in terms of a hypergeometric function was given in 
Ref. [5], and can also be found in the Appendix A. In this region the Coulomb interaction 
among the e'^e~ pair after soft photon radiation has a strength given by ~ and is 

resummed to all orders through the i ~ 1 Coulomb Green's function. Note that the usual 
QED perturbative expansion does not accommodate binding effects and cannot be applied 
to obtain the photon spectrum in this region. For smaller photon energies, u <^ ma^ the 
parameter v ^ 1 and tteiuj) — 2uj /ma^, in agreement with Low's theorem prediction for 
this decay [1]. 

• cj ~ ma {k, ~ \/rrujj): since now k ^ 1/a, the falloff of V{ijj) is ruled by the exponential 
factor from the Green's function, so the characteristic distance is r ~ k^^ ~ l/ma^/^. 
The argument of the trigonometric functions is constrained to ujr ~ \/uj/m ~ a^/^, which 
allows to multipole expand the full result also for radiated photons of this energy. The 
photon wavefunctions can be expanded out from the beginning and the leading term in 
such expansion is again given by de{oj) above, although higher terms yield now corrections 
proportional to {cu/m)"': 

= 4M (1 + 0(a;/m)) . (22) 

A new feature arises in the whole region ma^ <^ u <^ m because the parameter i/ ~ f 
is small and allows to further expand the Coulomb Green's function inside de{u!): 

G,{0, r-K) ^ ^ e- ^ G({0, r; k) , (23) 

i.e. the leading term of the expansion in u of the P-wave Coulomb Green's function is equal 
to the P-wave projection of the free Green's function. The latter shows that we can treat 
the Coulomb interaction as a perturbation in this region, and higher order terms in the 
i^-expansion of Gi correspond to insertions of the Coulomb potential in the diagrammatic 
picture. Using the result (23) and expanding out the o-Ps wave function e"''/", one obtains 
the leading term in u of de{uj), 

dJu) = -— / drr^G{iQ,r]K) + ... = --v + Oiv'^) 
ma J 3 

2a [m ( Lo mo? \ r,, , 

= -tVi:0-8;s-i^+-)^°<'''' '''' 

which agrees with the i/-expansion of the exact result of de{uj), which can be found in Eq. (A9) 
of the Appendix A. Kiojm correction to the lowest order result is generated automatically 
in the Oiv) term of Eq. (24) from the recoil correction which is included in k. Higher order 
terms in the v expansion correspond to the neglected terms in the expansion of G\ and the 
o-Ps wavefunction. The next term in the expansion of the photon wavefunctions, (ii^'*(a;) 



^ One has to approximate k ~ \/w — m£^o to make contact with the result in Refs. [1, 5], since recoil 
corrections were not considered there. 
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in Eq. (20), also produces a leading term proportional to v{uj/m) ~ a/2-\Juj /m, which is 
obtained by replacing Gi — > G{ and e"''/" 1: 



41) (^) ^ ^ / (irr3Gi(0,r;/.)e-'-/«^^''^ 
ma J 



2 



40 



5ma 



1 drr^ G{(0,r; «;) + ... = 1^ (^)V 0(^2^) . (25) 



The relation 



can be used in Eq. (25) to write d^e\i^) as an expansion muo/m. The exact result for d^e\i^) 
can be found in the Appendix A. For completeness we also give there the exact result of 
the full electric amplitudes 'D{ijj). 

Corrections of 0{(jj/m) can also arise considering higher powers in the momentum expan- 
sion of the annihilation amplitude A'^^^^ The higher order terms in p will pick intermediate 
e^e~ pairs with higher angular momentum L when they are transformed into derivatives 
acting on the Green's function. The scaling of this part of the electric amplitude is of order 

L-l / \ L-1 



(0|p(P) " G'|x)~(-) " (0|pG|x), 

Vm/ Vm/ 



i.e. shows a suppression ~ {u/m)^^ with respect the L = 1 case, since p ~ k ~ -Jujm. 
On the other hand, the radiative transition to the intermediate e^e~ L-wave state, ^5*1 — > 
7 + ^Lj, will select a L-tensor in the multipole expansion of the photon wavef unctions, 
that will yield a {ujr)^~^ ~ {ujjm)~^ extra factor with respect the lowest order transition 
amplitude [L = 1). Combining both factors we see that higher order terms in the expansion 
of ^(2'^') for the electric amplitude are thus down by an overall factor {ujjm)^~^. In order to 
compute these corrections it is convenient to write the expanded amplitude of A^'^'^'^ in terms 
of the spherical harmonics built from the vector p, so they can project the required L-wave 
component of the Green's function directly. A representation of the spherical harmonics in 
terms of the Cartesian coordinates of the vector is specially suited for this task, and can be 
found in Ref. [9]. 

Let us now turn to the magnetic contribution. We will not detail all the steps as we did 
for the electric amplitude, but just quote the result: 



M 



m 



-t 



n,P 

e 



(0| |n, P) i (n, P| i/x {cr^ ■ B(xi) + ■ B(x2)} |o-Ps) 

Eo — En,p — LU 



with 



-i- V'o(O) Wo 6 ■ (xV Aiu) , (27) 
m 

A{u;) J drr Go{0, r; k) e"''/" sin (^^) , (28) 

and 5 = (k X e)/a;.. Photon emission through a cr ■ B term can induce a transition from a 
triplet to a singlet spin configuration. The constant term in the annihilation amplitude A^'^'^'> 
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allows only for intermediate e^e~ states with quantum numbers ^5*0, which are contained in 
the £ = component of the Green's function. As in the case of the electric amplitude, the 
falloff of the Green's function fixes the characteristic distance for the multipole expansion to 
r ~ K"^, so that we can always expand in ujr for cu <^ m. Under the dipole approximation, 
the photon wavefunctions are dropped and the radiative matrix element reduces to the 
scalar product between the ground state and the rest of Coulomb wavefunctions, which are 
mutually orthogonal. In that limit, only the ^5*0 (p-Ps) ground state contributes [1], whose 
energy differs from that of the ground state o-Ps only if the hyperfine splitting is considered. 
The latter is only relevant for the study of the o-Ps decay spectrum at very low photon 
energies, u ~ ma"^ [1], and can be neglected here. The multipole expansion of the magnetic 
amplitude A{u!) then reads: 



-r/a 



(curY 



+ ... 



A{(jj) = 47ra; J dr Go{0,r; k) . ^ 

= am{uj) + a!;^\u) + . . . 

= am(u^)(l + 0{u;/m)) , (29) 
where the magnetic amplitude in the dipole approximation, 

1 + uj/Am 

is trivially equal to one if recoil corrections are neglected. If we take into account higher 
terms in the expansion of the photon wavefunction, transitions between the ground state 
and ^-wave radial excitations are possible, and the Coulomb interaction becomes relevant. 
For photons with energies in the range ma^ <^ cu <^ m the Coulomb interaction can be 
treated perturbatively and the first 0{a) correction arises in the subleading term in the 
multipole expansion, am{uj), and is proportional to v{uj/m) ~ a/2^uj /m (see Eq. (A15)). 

The exact result for a"m (^) and for A{uj) are given in the Appendix A. 

According to Eqs. (20,29), the NRQED calculation of the o-Ps decay amplitude at leading 
order in the multipole expansion reads: 

M^Me + Mm^-e V'o(O) (^(x^ Wi • e 0)^0 a,{uj) + ^ ^ ^5 • (x^ (t>)e, a„.(a;)^ . (31) 

From a knowledge of the amplitude, we can derive an expression for the o-Ps differential 
photon spectrum at low energies (a; <C m) written in terms of the electric and magnetic 
amplitudes [1]: 



dr ma 



6 



dx 97r 



-X 



7 



|2 I _|_ I |2 



X = uj/m , (32) 



Let us now focus on the region oo ~ ma. It is clear that since in this region Coulomb effects 
can be accounted for perturbatively, one can apply the conventional QED loop expansion. 
The leading (LO) and next-to-leading (NLO) terms in a of the o-Ps 37 NRQED amplitude 
when X ~ a are obtained by taking 

2 a 

ae = l--^ + 0(Qi) , 0^ = 1 + 0(0;) , (33) 
o yx 
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and the NRQED spectrum in this region, up to NLO, thus reads: 



dr 

dx 



2ma^ 




(34) 



277r 



The LO term of the NRQED amphtude (31) matches the uj ^ hmit of the tree-level 
QED calculation first done by Ore and Powell [10], as it was shown in Ref. [1]. The same 
is of course also true for the LO photon spectrum. The comparison for the a^/x correction 
requires the one-loop result of the QED series. In Sec. Ill we will derive the soft photon 
limit of the 1-loop QED spectrum obtained from an analytic expression of the phase-space 
distribution given in Ref. [6], and check that the leading term when a; — > agrees with the 
ay/x correction in the NRQED spectrum shown above ^. In Sec. IV we will compute the 
leading term of the 1-loop QED amplitude for soft photon energy by means of the threshold 
expansion, and show that it arises from a loop-momentum region p ~ y/um, and that it 
indeed agrees with the NRQED result. Concerning higher orders in the uu/m expansion for 
the 0{a^) and 0{a) terms discussed above, any possible mismatch between the QED and 
the NRQED amplitudes would require that we introduce a short-distance e'^e" 37 local 
contribution in the NRQED Hamiltonian ^. 

III. ONE-LOOP QED SPECTRUM FOR ma^ < a; < m 

An analytic evaluation of the 1-loop o-Ps — > 37 decay amplitude has not become available 
until recently [6]. In the later work, Adkins has provided a compact form for the one- 
loop phase-space distribution, which allows to obtain the one-loop correction to the energy 
spectrum upon integration: 



where Xj — uji/m{i = 1,2,3), with ujj the energy of the photons, which satisfy that Xi + 
X2 + X3 — 2. Analytic expressions for the function F{xi,X3) can be found in the Appendix 
of Ref. [6]. They involve a number of intermediate functions of the variables Xi, including 
dilogarithm functions. Therefore an analytic integration of the formula (35) may be difficult 
to obtain. However, for an evaluation of the low-energy spectrum one does not need to 
perform the full integration before taking the limit — > 0. A Taylor series in xi can be 
obtained by first taking the small photon energy limit of the integrand in Eq. (35) up to 
the desired accuracy, and afterwards doing the (trivial) integration. Since also the variables 
I — X2 and 1 — ,T3 are small when Xi ^ , some care is needed in order to perform this limit. 

The two essential variables that parameterize the 1 — > 3 phase-space have been choosen 
in Eq. (35) to be xi and X2, where xi is the variable of the differential spectrum. Therefore 
xi <^ X2 1- Note that in the X2 integration the limits imply that 1 — X2 < xi, so we can 



* A comparison between the effective theory 0{a) spectrum and a numerical calculation of the one-loop 

QED result [11] has been done by Voloshin [5], showing a good agreement. 
^ Let us remind that the NRQED Hamiltonian is constructed to reproduce the QED amplitude neglecting 
binding effects. 




(35) 
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write 1 — X2 = xit, with t e [0, 1], and change the integration variable to t: 



•ir- - f\t-, pi^l^ ^, (36) 



dxi Jo {l-xit){l-xi{l-t))' 

where 

F{xi, t) — F{xi, 1 — Xit) + permutations . (37) 

The integrand in Eq. (36) is symmetric under the change t ^ {1 — t), which is indeed 
equivalent to changing X2 <-> x^- With the redefinition above we can Taylor expand the 
integrand in the variable Xi while keeping the variable t fixed, and the resulting terms turn 
out to be well behaved in the Xi limit. Up to the second non-trivial order it yields a 
relatively simple phase-space distribution: 

F(xi,t) IGtt , , 

^ ^ - [3-4t(l-t)]^ 



{1 - xit){l - xi{l - t)) 3 

+ [2t{l-t) (37r^ - 48 log 2 + 40) + tt^ + 16 log 2 - 60 ] xi 

+ 0{xf), (38) 

that can also be written back in terms of X2 if needed. Accordingly, the 1-loop QED 
correction to the low-energy spectrum reads (writing Xi x) 



dFi 2ma^ 



dx 277r 



14 ^ a /37r^ 35 \ o/o, 
-—a^ +-l — - — )x + Oiax'/^) 



(39) 



We see that the a^/x term above agrees with the corresponding one in the NRQED formula, 
Eq. (34). The ax correction also matches the corresponding term in the EFT computation, 
which was included in the formula (34) by Voloshin [5] considering the 0{a) correction 
to the electron gyromagnetic ratio and the 0{a) corrections to the annihilation ampli- 
tude e+e~ — > 27. The latter correction can be accounted for in a straightforward way if 
the NRQED spectrum is written in terms of local four-fermion operators which effectively 
integrate out all possible states of hard photons [1]. In this description, the coefficients mul- 
tiplying the squares of the electric and magnetic amplitudes in the photon spectrum (32), 
are identified with the imaginary part of the matching coefficients of S- and P-wave four- 
fermion annihilation operators (see Eqs. (21,74) of Ref. [1]). These coefficients can be found 
up to 0{a^) in Ref. [12], and yield a part of the ax correction shown above. The remaining 
part, due to the radiative correction to the electron gyromagnetic ratio, is accounted for in 
the NRQED Lagrangian naturally as the 1-loop matching coefficient of the a ■ B interaction. 

To conclude, we wish to point out that higher order corrections in the x-expansion of the 
1-loop low-energy photon spectrum can be computed straightforwardly from the analytic 
phase-space distribution given by Adkins through the expansion procedure outlined here. 



IV. THRESHOLD EXPANSION OF THE QED 1-LOOP DIAGRAMS 

The threshold expansion technique is a prescription developed by Beneke and Smirnov [7] 
for the asymptotic expansion of loop integrals involving heavy massive particles close to 
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FIG. 2: Ladder and double- vertex graphs contributing to the o-Ps decay amphtude to order a. 
Momenta flow in the same direction as the arrows. 

threshold [i.e. that are moving at small velocities). It provides the velocity expansion of an 
amplitude up to a certain order as a set of simpler integrals than those present in the original 
Feynman diagram. The method relies on the so-called strategy of regions, which replaces the 
integration over the whole space of loop momenta to the integration only over some specific 
regions. The first and crucial step is to identify the relevant momentum regions in the loop 
integrals. These in principle follow from the pole structure of the Feynman propagators, 
that is characterized by the relevant scales that appear in the problem. For amplitudes 
involving a non-relativistic pair, three scales enter the dynamics: the mass of the particle, 
m, the relative 3-momentum ~ mv and the non-relativistic energy ~ mv^ . Accordingly, four 
different loop momentum regions were identified to give non-zero contributions in Ref. [7]: 

hard : go ~ q ~ ? 

soft : go ~ q ~ '^'v , 

potential : go ~ q^/m mv'^ , 

ultrasoft : go ~ q ~ mt>^ . (40) 

In a second step the original integration region must be decomposed into a sum of inte- 
grals, one for every region above, and a Taylor expansion in the parameters that are small 
in each regime performed. Every integral, containing just one scale, thus contributes only to 
a single power in the velocity expansion, which can be determined easily before integration. 
The procedure requires the use of dimensional regularization, even if the original integration 
is finite ^. Following these heuristic rules, the authors of Ref. [7] reproduced the exact v 
expansion of some one-loop and two-loop examples. Although a formal proof of the validity 
of the asymptotic expansion close to threshold has not been given, the perfect agreement in 
the examples supports their use in general one-loop diagrams. 

We now wish to apply the threshold expansion to the 1-loop o-Ps 87 amplitude 
when one of the photons has an energy u ~ ma, where a conventional perturbative QED 
calculation is still valid. The momentum of the emitted soft photon, k ~ ma, introduces a 



^ Other analytic regularization procedures can be employed, since the important point is that scaleless 
integrals are put to zero. 
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further scale in the amphtude and, as we will show below, gives rise to a new region that 
contributes to the exact result. It corresponds to a loop momentum 



(41) 



which arises from the heavy particle propagators and indeed resembles the potential region 
but with a different scaling for the 3-momentum q ~ ^Jmio. It thus corresponds to a small 
momentum region {q <^'m) which lies in between the hard relativistic region and the small 
momentum regions shown in Eq. (40), and we shall refer to it in the following as 'soft- 
radiation' region. The physical origin of the new region has been discussed in the NRQED 
description of the o-Ps decay in Sec. II. 

The 1-loop graphs contributing to the o-Ps — > 87 amplitude can be classified into one 
of the following categories [6, 13]: self-energy, vertex corrections, double- vertex, ladder and 
annihilation contributions. The leading term in the uj/m expansion is given by the ladder 
and double-vertex graphs, shown in Figs. 2a and 2b, respectively. The momenta flowing 
through the fermion propagators have been written in the figure. The Positronium atom 
four- momentum is P — ki + k2 + k^, with components P — (M, 0) in the o-Ps cm. frame, 
and we can take M ~ 2m safely for the discussion that follows. Note that the assignment 
of momenta in the figure implies that the electron and positron inside the bound state 
have incoming momenta P/2 = (m, 0). The evaluation of the 1-loop diagrams for the o-Ps 
decay with e'^e" strictly at threshold is a valid procedure to obtain the matching coefficients 
of the NRQED annihilation operators that give the o-Ps total decay rate, as it has been 
shown in Ref. [13]. Setting the velocity to zero automatically kills off the contributions of the 
small momentum regions shown in Eq. (40), so that only hard and soft-radiation momentum 
regions will survive in the QED 1-loop o-Ps amphtude ^. In order to show that the threshold 
expansion technique for diagrams with a soft photon emission does not only apply for this 
particular kinematic configuration, we will consider in the Appendix B a more general case 
with P^ 7^ 4m^. where several small momentum regions exist at the same time. 

Let us consider the ladder graph shown in Fig. 2a, where the soft photon is the one 
attached to the incoming electron. Therefore ki ~ ma, while k2, k^ are hard photons 
scaling as m. Two integrals turn out to be relevant in the a;i — > limit of the ladder 
amplitude. The first one is the 5-point scalar integral: 



?2 [(q + P/2f - m2] [{q - P/2)2 - m^j [{q + P/2 - k^f - m^j [{q - P/2 + k^f - m^j ' 



The integration measure is defined as [d^q\ = d^q/i{2Tr)^, m. D — A — 2e dimensions, and 
the standard -\-ie prescriptions are implicitly understood in the propagators. When the loop 
momentum is hard, we can expand the fermion propagator in ki. In this way, only the scale 
m enters into the propagators in the hard region. Taking into account that [d^q\ ~ m^, 
the leading term in the contribution from the hard region is estimated to be /q'^'* ~ m~^. 
However, for tui — > the dominant contribution to the integral above comes from a loop- 
momentum region where q <^ m. Since the momentum that fiows through the fermion 



Loosely speaking, in the approach of Ref. [13] the non-relativistic dynamics is properly taken into account 
by the NRQED part of the calculation. 




[d^q] 



(42) 
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propagator connecting the two hard photons is always hard, when q <^ m we can expand 



''■^-^'■*^. + ... (43) 



[(5 - P/2 + ksf - m'] -P -k, (P ■ k,)' 

so that the propagator effectively shrinks to a local annihilation vertex. Likewise, in the 
remaining massive propagators we can drop the terms with respect q ■ P — 2mgo- The 
leading contribution of /q in the small momentum region is thus 

/(^-'^u) ^ L_ f . (44) 

2mus J {ql — q^) (2mgo — q^) {~2mqo — q^) {2mqo — 2mu!i — q^) 

We have also dropped the term q ■ ki <^ 2mgo- Closing the integration contour in the 
upper complex qQ-plane, we pick the contribution from the pole of the massive propagator 
at go = — q^/2m + ie. The contribution of the residue of this pole to j^^™^^^^ reads 

with n = D — 1. In deriving Eq. (45) we have also expanded the photon propagator in q^, 
since <^ q^ once we pick the residue from the fermion pole. The integral in Eq. (45) is 
finite and it is dominated by loop momentum q ~ ^/mlJ, which corresponds to the scaling 
of the soft-radiation loop momentum region introduced above ^. As in the hard region, 
the contribution from the soft-radiation region to the final result can be readily estimated 
before making any integration by first expanding the integrand according to the hierarchy 
go ~ q^/m ~ muJi. For the leading term the propagators yield a factor ~ m~'^{muJi)~^, 
while [d^q] ~ m~^(q^)^/^ ~ m"^(ma;i)^/^, so that we obtain an overall scaling ~ 
m~^(mcc;i)^^/^, as found in Eq. (45). The exact result of the integral Iq can be easily derived 
from the 4- and 5-point functions computed in Ref. [6]. We have checked that the leading 

(s-r) 

term in the a;i — > limit of the exact expression for Iq agrees with 7q above. There is 
also a contribution to j^^™^") from the pole of the massless propagator at go = — |q| + ^e. 
After picking up the corresponding residue the remaining integration is dominated by loop- 
momenta |q| ^ LJi. The characteristic loop- momentum in this region is thus the same as 
for the soft region in Eq. (40), if we consider v a. The contribution from this pole has 
an overall scaling [d^g]/(mgo)^ ~ m~^6i;f^, i.e. larger than the hard region one, although 
suppressed by ^JmLJl with respect the soft-radiation region. We do not give the result 
because we will not keep terms of the similar order that arise from other diagrams. 

The second integral that we shall need for the leading term in the a;i — > limit of the 
ladder amplitude is the 4-point tensor integral 

^ I \jD -\ / .px 

° i (f [(g + ^/2)2 - m2] [(g - P/2)2 - [(g + P/2 - kxf - m^] ' ^ > 

with the indices i, j = 1, 2, 3. We have already expanded out the hard fermion propagator in 
order to write down /q-' , since the leading contribution will come up from the g -C m regime. 
As we shall see when we put back the numerators of the amplitudes, part of the g'g-' term 



The integral in Eq.(45) with wi replaced by Aqcd was also found in a non-perturbative enviroment [14]. 
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arises when a is brought to the numerator of the integrand from the subleading term in 
the q <^m expansion of the hard propagator shown in Eq. (43). Strictly, one should do the 
tensor reduction of the integral and afterwards apply the threshold expansion to the scalar 
integrals that appear. However, since the tensor structure in the integral in Eq. (46) only 
involves the spatial indices, we shall take the path of expanding the integrand according to 
the different loop-momenta directly without prior reduction to scalar integrals. In the small 
loop-momentum region, the propagators in identical to those in j^^™^")^ thus 

they share the same pole structure. The expansion of the propagators in the soft-radiation 
region leads to: 

^^•(s-r) = J_f un ] y ^ = -J: ^ + 0(e) . (47) 

After dropping the ki-dependcncc in the fermion propagator, the reduction of the tensor 
integral is trivial. The contribution from the gluon pole in Iq can be shown to scale as 
~ cui/m^, and it is thus suppressed with respect the soft-radiation contribution. 

For the double vertex diagram of Fig. 2b the leading term in the a;i ^ limit comes 
from the integral 

J ^ f .4o^ 

J q''[{q-P/2y-m^][{q + P/2-kiy-m^]' ^ ^ 

after expanding the hard region propagator as done in Eq. (43). The expansion of Ii in the 
small loop-momentum gives an analogous expression to (44) but with one fermion propagator 
less. The soft-radiation region arises from the fermion pole at qo — — q^/2m -|-ie, and yields 

.(s-r) i /" r •, 1 1 , ^_l 



4m J (q"' + mwi) lonm 

On the other hand, the contribution from the gluon pole in the small loop-momentum scales 

as m~^, so it is again suppressed with respect the soft-radiation region. 

Let us now include the numerators of the propagators and vertex factors in order to 
calculate the full amplitude. We can use the identity 

^ + m _ Es^s(p)M.(p) , J2s^s{-P)vsi-P) 



p2 _ jj^2 _|_ pO _ _^ pO _|_ _ 



with Ep = \/ + p2 ~ m + p^/(2m), to pick the numerators associated to the non- 
relativistic fermion poles that arc relevant when the loop-momentum is small, and after- 
wards perform the non-relativistic expansion of the spinors in p/m. This will automatically 
lead us to the usual non-relativistic Feynman rules for the vertices and propagators. The 
identity (50) effectively separates particle and antiparticle contributions, as it is conven- 
tional in the time-ordered perturbation theory. The spinors in Eq. (50) have non-relativistic 
normalization, more precisely: 



where s is the spin label. Let us illustrate how the non-relativistic expansion is performed for 
the ladder diagram. The structure of the fermion propagators in the small loop-momentum 
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M,;(q - ki) 

tv(-q) 



+ perm. 



FIG. 3: The soft photon vertex (A^i) and the 2j annihilation (A^23) parts in the ladder graph. 

region has been shown in Eq. (44). Take, for example, the fermion propagator at the left 
of the soft photon vertex, with momentum q + P/2. The expansion of this propagator for 
q <^ m yields 2mg° — q^. The latter agrees (up to a normalization factor) with the first de- 
nominator that is obtained from the identity (50), which reads q^ + /2 — ~ g° — q^/2m. 
Therefore we will associate the electron spinors Ws(q)'Us(q) to this denominator in order 
to derive the leading term in the small momentum region. In the same way, the expanded 
fermion propagators after soft photon emission and after the vertex correspond to electron 
and positron contributions, respectively, in the splitting of Eq. (50). For the massive line 
connecting the /c2 to the vertices, we need to keep the full covariant propagator because 
the momentum flowing through it is always hard. Taking into account the particle and 
antiparticle spinor wavefunctions from the adjacent propagators, the non-relativistic reduc- 
tion of this part of the ladder diagram is indeed equivalent to the e~^e~ — > 27 annihilation 
amplitude of an incoming electron with 3-momentum q — ki and an incoming positron with 
3- momenta — q (see Fig. 3b), 

M23 = -I Vs\-ci) /a _^ i2 M.(q - ki) + {2 ^ 3} , (52) 

where we have also included the contribution from the graph with the hard photons per- 
muted. Recall that we found already the non-relativistic expansion of the e'^e" — > 27 
amplitude in Sec. II. Indeed, the expansion of Eq. (52) for the hierarchy m ^ q ^ ki, char- 
acteristic of the soft-radiation region, yields the same leading and subleading terms already 
shown in Eq. (5), so that 

-M23 = xl' (Wo + Wi ■ q) + 0{a) . (53) 

The size of the terms in the expansion of M.2z above should be estimated considering the 
soft-radiation region scaling, q ~ ^JmUl ~ ma^/'^. The reason why we only need to keep 
the two first terms in A^23 shall become clear when we put together all the pieces in the 
amplitudes. Similarly, the non-relativistic expansion of the soft-photon vertex is also done 
easily once we pick the particle contribution from the propagators (see Fig. 3a): 

A^i = -i e M,(q - ki) ^iM^(q) = q ■ ei 0^0^ + — (ki x ei) 0j<T0^ + C(a^/2) . (54) 

m 2m 

We can check that the non-relativistic reduction of the amplitude M.i agrees with the 
Feynman rule that is derived for the e~ — e~ +7 vertex from the non-relativistic interaction 
Hamiltonian, Eq. (1). Finally, the photon exchange between the initial e^e~ reduces to a 
Coulomb-like potential: 

= —Ur{qh^u{0)v{0)l,vA-cD = --^(PUx^Xs' + 0{a') , (55) 
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in the region go jva ^ m, and with 0. x the Paiili spinors of the initial e^e~ . We 
can now put together the different pieces, M.^, M.^ and 7^23, including the denominators 
and taking into account that ^s4>s4>\ — SsX^Xs — 1- The leading term of the threshold 
expansion of the ladder amplitude in Fig. 2a (plus the graph with the photons ki and k-^ 
interchanged) is thus equal to 

(x^Wi-q^) (q.e)-2Wo^i-(xto-0) 
(q^) (9o - G^jlm) {qo + qV2m) (go - cui - qV2m) 

= e x^Wi-e(/. -W— + ^e^5i- X^fr0)-./— . 56 

The integrals appearing in Eq. (56) are proportional to the soft-radiation contributions of 
the integrals Jo and Iq , computed in Eqs. (45) and (47), respectively. We can proceed in 
a similar way to obtain the leading term in the threshold expansion of the double-vertex 
diagram, Fig. 2b: 

2m J (q^) (go + qV2m) (go - Ui - q^/2m) 

= -ie^«..(;,V0)"/^. (57) 

m 2 \ uji 

where we have required the contribution from the soft-radiation region of /i, Eq. (49). For 
the double- vertex amplitude there is no contribution from the q ■ ei term in the soft-photon 
vertex because the incoming electron is static. The Wo terms fr om th e ladder and double- 
vertex graphs cancel each other and we are only left with a a^mjujx proportional to Wi. 
It is straightforward to verify that the contribution from the ladder graphs where the soft 
photon vertex is attached in the lower part of the diagram {i.e. to the incoming positron 
line) gives the same Wi term. Therefore, the leading term in a;i ^ limit of the 1-loop 
e+e~ — > 87 QED amplitude reads 

A^i-ioop = e(x^W,-60)^./^, (58) 

which agrees with the leading Oia) term in the NRQED computation of this amphtude, 
Eqs. (31) and (33), once we include the o-Ps wavefunction factor. 

A comment on the rest of 1-loop diagrams that had been not considered is appropriate. 
First one should note that all diagrams where the soft photon k\ is placed among the two 
hard ones scale as 0(1) because the momenta flowing through the fermion propagators that 
connect two photon vertices are always hard. With respect the annihilation contribution, 
the scaling can be inferred from that of the ladder graph if we replace the t-channel photon 
exchanged between the initial e^e~ by a s-channel photon exchange. The replacement 
carries a q^/m^ suppression in the non-relativistic expansion, which is of order uji/m for the 
soft-radiation contribution. Regarding the self-energy and vertex correction topologies, they 
correspond to typical relativistic corrections whose effects are not expected to be enhanced 
when we explore the non-relativistic regime. Indeed it can be easily checked from an quick 
analysis of the propagators in the small loop-momentum region that for the integrals in the 
vertex and self-energy corrections, we can always avoid the fermion poles so that only the 
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(a) 






27 



FIG. 4: (a) NRQED diagrams that contribute to the 0{ay^m/uJi) term of the e~^e~ — > 87 decay 
amphtude. The zigzag hnes represent soft photons; {b) local vertex that describes real soft photon 
radiation from the heavy fermion-antifermion pair in vNRQED. The 27 label denotes that the 
vertex arises from the integration of two hard photons. 

pole in the photon propagator at g° ~ q contributes. In this way, the small loop-momentum 
expansion of the 2- and 3-point scalar integrals give results which are down by {ui/m)'^ and 
uji/m, respectively, with respect the hard region contributions to the same integrations. 

It is interesting to note that the expressions (56,57) correspond to the amplitudes that 
would be written down from the diagrams displayed in Fig. 4a using the usual NRQED 
Feynman rules with free fermions as asymptotic states ^, and after we drop the ki term in 
the fermion propagator: 



Keeping the full NRQED propagator and expanding afterwards leads of course to the same 
leading order term, but unnecessarily complicates the calculation. In the usual NRQED 
counting for bound state systems, the heavy fermion and antifermion inside the bound 
state have energies E ~ ma^ and relative 3-momentum q ~ ma, so they are classified as 
potential degrees of freedom. As we explained in Sec. II, the radiation of a soft photon from 
the e~^e~ allows for higher virtualities in the relative 3-momentum, of order q ~ y/muji. 
Therefore the heavy particles propagating in the loops of Fig. 4a do not obey the potential 
scaling (otherwise the massive propagator (59) for a potential fermion will reduce to just 
l/ui and we will miss the contribution from the soft-radiation region). The fact that we 
have to further expand the NRQED propagators in the soft-radiation regime is an unwanted 
feature because it jeopardizes power-counting. However, we can take advantage from the 
fact that soft radiation takes place at distances ~ 1/ ^/muJ, smaller that the typical size of the 
Coulomb bound state, to integrate out the contributions stemming from the soft-radiation 
momentum region. From the point of view of the potential fermions, soft radiation is a 
short- distance process, just like the 27-annihilation, and can be described by a local vertex. 
In the vNRQED/vNRQCD [15] formalism (with the external sources for the hard photons 
integrated out) the decay spectrum in the soft-energy can be calculated from the imaginary 
part of the matrix element of a 6-field operator as the one shown in Fig. 4b. The matching 



^ Recall that we used bound state perturbation theory in coordinate space to compute the NRQED am- 
plitude in Sec. II. In the w ma region the leading 0(a) correction in the latter computation could be 
obtained by replacing the intermediate Coulomb Green's function by the free one and neglecting the o-Ps 
wavefunction dependence (see Eq. (24)). This is effectively equivalent to consider the Coulomb interac- 
tion as a perturbation, so that in this regime the NRQED computation can be done using as unperturbed 
states either free propagating fermions (as in Fig. 4a) or e+e^ Coulomb eigenstates (as in Sec. II). 



1 1 



(59) 



g° — tUi — (q — ki)2/2m q'^ — Ui — /2m 
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coefficient of such operator can be obtained by comparing with the QED result in the 
uji ~ ma region: at tree- level it is just a constant xa^ (the limit a; ^ of the Ore-Powell 
decay spectrum [10], divided by the x factor from the phase-space), while at 0{a) it is 
proportional to a^^fmfu). Note that this picture is similar to the NRQCD factorization 
approach [12] used to describe the production and annihilation of heavy quarkonium, that 
we briefly discuss in the next section. 



V. RELEVANCE FOR HEAVY QUARKONIUM 

The threshold expansion technique applied above for the o-Ps— > 87 decay, can also be use- 
ful to calculate the direct photon spectrum of quarkonium radiative decays in the soft-energy 
region. Heavy quarkonium decay rates can be understood within the NRQCD factorization 
approach [12], which allows to separate the short distance physics related with the heavy 
quark-antiquark annihilation process from the long distance bound state dynamics. An 
operator product expansion can be written down for the direct photon spectrum in, for 
example, T decays: 

Ay 

- = 5]C„(x)(T|a|T), (60) 

n 

where the sum above extends over all (5Q['^] configurations that can be found inside the 
quarkonium, and the C„(x) are short distance Wilson coefficients that can be determined 
from the annihilation cross section of the on-shell QQ\p\ pair as a perturbative series in 
^^(mfe). The C„ in the long-distance matrix element are NRQCD operators which are 
organized in powers of the relative velocity of the heavy quarks. At leading order only 
the color-singlet operator 0\(^S\)^ that creates and annihilates a quark-antiquark pair in a 
color-singlet '^S\ configuration, contributes. The nonperturbative NRQCD matrix element 
is related to the T wave function at the origin 

(T|Oi(35i)|T) = (T|(/>t,o-VpxLp,tTVp'|T) = 2A^,|V.(0)|^ (61) 

where now 0p, x_p refer to heavy quark and antiquark fields respectively. The matching 
coefficient at leading order in ctg comes from the tree-level QQ — > gg^ annihilation [16], 

2 — x (\ — x)x 2 (1 — x)^ log(l — x) 2 (1 — x) log(l — x) 
X ^ {2-xf (2 - xf ^ J ' 

(62) 

which has identical x-dependence to that of the Ore- Powell spectrum for o-Ps [10]. In- 
deed the factorization approach applies trivially to QED bound states, where everything is 
calculable perturbatively and the bound state wavefunction reduces to the Coulomb one. 
The leading order coefficient for Ps decays, C\[x\ is obtained by the trivial replacements 
alQl/ml o? jn? and ?>2/21 — > 16/9 in Eq. (62). In the soft-energy region, it reads 



Ci(x) 



32 aalQl 
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The 0{a) correction to the o-Ps— > 87 amphtude for soft photon energies computed in the 
previous sections can also be included, so that 



14 

5x — — a^/x + 0{ax^) 



(64) 



that we can plug in Eq. (60) to obtain the o-Ps photon spectrum to NLO in the soft- 
energy region (we also need to replace — V'c and Nc — > 1 in Eq. (61)). A similar result 
to that of Eq. (64) could be obtained for the 0{aa^) short-distance coefficients C„(a;) in 
heavy quarkonium decays by doing the threshold expansion of the corresponding 1-loop 
QCD diagrams. The latter could provide a determination of the 0{as) corrections to the 
quarkonium photon spectrum in the soft-energy range of photon energies. To our knowledge, 
as corrections to the coefficient Ci{x) in Eq. (62) are only known numerically [17]. We 
should be aware though that the fragmentation contributions to the photon spectrum in 
quarkonium decays {i.e. those in which the photon is emitted from the decay products) 
become important in the \ow-x region [18, 19], namely for x < 0.3 in T decays, and have to 
be properly taken into account for a computation of the full spectrum. At lower energies, 
X < a'^, the emission of the photon can produce transitions to virtual bound states and 
does not longer belong to the short-distance part of the decay [1]. The standard NRQCD 
factorization also breaks down at large values of the photon energy [20], where one needs to 
consider also coUinear degrees of freedom [21]. 



VI. CONCLUSIONS 

The region of photon energies uj ~ ma in the radiative decay of a heavy fermion- 
antifermion pair has been analyzed in this work through the study of the thrcc-photon 
decay of ground state orthopositronium. The soft-energy region effectively separates the 
regime where binding effects become essential, from the hard-energy region, where the de- 
tails of the bound state dynamics are irrelevant. The NRQED framework is able to yield 
the correct o-Ps spectrum from energies u < ma"^, where it properly accounts for bound 
state effects, up to energies co ~ ma, where the binding can be neglected and the calculation 
is simplified by replacing the Coulomb Green's function with the free one. We have found 
agreement in the comparison of the o-Ps spectrum in the soft-energy region as computed 
with NRQED and with conventional perturbation theory in QED. For the NRQED compu- 
tation, the leading order approximation in the multipole expansion of the photon field yields 
already the dominant contribution in the whole energy range a; -C m, as it had been argued 
by Voloshin [5] . We have explicitly calculated the contribution from higher order multipoles 
to the o-Ps decay amphtude in NRQED, and discussed the size of relativistic corrections 
coming from other sources in the soft-energy region. 

The 1-loop QED result of the o-Ps spectrum in the soft-energy region has been obtained 
by two different methods. First by expanding the analytic expression for the phase-space 
distribution obtained by Adkins [6], and second by means of the threshold expansion tech- 
nique. With the latter method it has been shown that the dominant contribution arises from 
a new loop-momentum region ~ q^/m ~ ma. This momentum region arises naturally 
when the heavy fermion-antifermion system decays emitting a soft photon and has to be 
added to the rest of momentum regions which are known to be relevant for the analysis 
of heavy particle- antiparticle loop diagrams. The identification of momentum regions that 
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have not been previously considered is mandatory for the success of the asymptotic expan- 
sion method in a wider range of kinematic situations. It is also of conceptual importance 
for the construction of EFT's for these systems with well-defined power-counting rules and 
able to reproduce the correct low-energy behaviour. 
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APPENDIX A: CLOSED FORMULAS FOR THE NRQED O-PS^ 87 AMPLI- 
TUDE 

We give first some useful formulae for the computations that follow. The partial waves 
of the Coulomb Green's function can be written in terms of the associated Laguerre poly- 
nomials [23] . Setting one of the arguments to zero they read, 



n=0 ^ / \ / 



r=0 

A related integral is 



that we shall need only for the £ = and I — 1 cases: 



a 



(Al) 



with [24] 

— n ■ V / 



JuM) = I dze-^Ll'^\z) = (-1)™^ , (A3) 



= -^--;^ + aia-n-3) + - + -n + 3. (A4) 

Let us start with the electric amplitude I^(a;), Eq. (14). Using the representation (Al) and 
making the change of variables z — 2Kr, it can be written as: 

viu) = 16. Q E;:t2^H^^"^ + £^-^^"^} + ^-'-' ^^^^ 



n=0 
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with 

1 + 1/ iu> 

a — . 

2 Ak 

(recall from Sec. II that u — l/na). With the help of the results (A3,A4) and the identity, 

E°° a:" ^ 2-^i(a, l;a + 
n+a~ a ' ^ ^ 

n=0 

we arrive to a closed expression for Vluj) in terms of a hypergeometric function, 

,//t\3r ^ Jq;(q;- 1) - :^^(2Q; + i/- 1) u ^ 
V{uj) = (-] Lf^{2 - j., 1; 3 - i/; /x)^ ^-i^^ ^ + 4^^" / + 

(A7) 

where we have defined /i = {a—l)/a. The expansionin {uj/kY of the result (A7) is equivalent 
to the one which is obtained by directly expanding the trigonometric functions inside the 
integral representation of Viuj), see Eq. (19). The leading term (dipole approximation) 
yields a representation of deioj), Eq. (21), in terms of an hypergeometric function, 

, 4z/(l + 2z/) 8 z/2(z/-l) ^, ^ ^ Z/-1 

The representation above agrees with the corresponding one for ae{uj) = 1 + de{uj) found 
in [5]. The expansion of de{u}) for small u is legitimate in the ma^ <^ u) <^ m region, and 
follows immediately from (A8), 

dM = + ^ (1 - log 2) + 0{v^) , (A9) 

where higher order terms in v account for further Coulomb interactions. The subleading 
term in the multipole expansion of the electric amplitude corresponds to the (uj/k)'^ term 
in the expansion of 'D{uj) and so on. Since {u/kY ~ uo/m (see Eq. (26)), multipoles of 
r-th order get suppressed by the overall factor {uj/mY . A formula for the r-th multipole is 
obtained in a more direct way through the a;-expansion of the trigonometric functions inside 
T>{ijj) prior to integration. It reads 



d^:\u) = (-l)'-+i-f±i- (^) 'y ^ +^(/3) , (AlO) 

' ^ ' (2r + 3)!\4«;/ ^ n + 2 - 1/ "'^ ' ^ ' 

^ ' n=0 

with = (1 + z/)/2. A closed expression for (ie\^') can be obtained with the aid of 
Eqs. (A3,A4) and a generalization of the identity (A6), 

oo r 

V ^^x" = V(-a)^-i (Li,_,(a;) + 5,,,) + (-a)'-^ 2^1(0, 1; a + 1; x) . (All) 

^^+« 

For example, the next-to-leading multipole (r = 1) reads 

,m/ X 1 (^\^ ^ rz/2-3z/-2 4z/(l + z/2) ^ ^ ^ 

2 



-5' 



Q +0{uW/^), (A12) 
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whose leading order term in the small v expansion agrees with our previous calculation using 
the free Green's function in Eq. (25). 

Let us also give the exact formula for the full magnetic amplitude ^(a;). From Eq. (28), 
we obtain 



oo ^ 

•4M = -|^E;rTT^<o(«) + h.c. 

n=0 



im 



1 + — 2^1(1-^^,1; 2 - z/; /x) ^ +h.c.. (A13) 



laK \ a{\ + v) 

As in the case of the electric amplitude, the multipole expansion of A{lo) is realized through 
a series in {oj/kY. The exact result for the magnetic dipole term, a^(a;), has been already 
shown in Eq. (30). The next-to-leading term in the multipole expansion reads 



(1)/ X ^ l_mi^ /ujy 1 

12 ^2 \,^) (1 -i.2)3(l + ^) 

X js - + + 5z/2 - 8z/2 2i^i(l - z^, 1; 2 - i/; ^^^) ) . (A14) 
The expansion of am (<^) for small v gives 



(1 - 1/2)3(1 + J,) 



where we have made use of the relation between k/^ and cu given in Eq. (26). Finally, we give 
the general form of a the r-th order magnetic multipole: 



^ ' n=0 

which scales as {uj/mY 



APPENDIX B: THRESHOLD EXPANSION OF A TRIANGLE DIAGRAM 
WITH A SOFT RADIATED PHOTON 

In this appendix we show a further example of the application of the asymptotic expan- 
sion method near threshold to a 1-loop diagram with massive lines and a soft momentum 
component radiated off. The diagram is shown in Fig. 5. This example is similar to the 
1-loop ladder diagram considered in Sec. IV, but we shall not take here the total incoming 
momentum equal to Aw?, but instead define y = w? — P'^/i ~ {maY. The parameter y, 
which appears in the massive propagators, will give rise to a non-zero contribution from the 
potential region, which was not present in the ladder diagram evaluated at threshold. For 
simplicity, we have substituted the Coulomb-ladder interaction by an effective production 
vertex. We shall take the soft photon momentum k ~ ma, and keep terms in the asymptotic 
expansion up 0{a). 
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FIG. 5: Scalar triangle diagram. 



If we consider scalar propagators, the integral of Fig. 5 is given by 

[d^q] 



(Bl) 



{q^ + q ■ P — y) {q^ — q ■ P — y) {q^ + q ■ P — 2q ■ k — P ■ k — y) 

where we have choosen a routing of the external momentum P through the massive lines of 
the graph which allows to use the scaling arguments defined in Eqs. (40) and (41). When 
the loop momentum is hard (g ~ m), we can expand the propagators in the small variables 
y and k. Keeping terms up to 0{a) the expansion in the hard region reads: 

.(h) ^ f Kg] A , 2q.k + P.k 

' J {q^ + q-PY{q^-q-P)\^ q^ + q ■ P ^' 

^ ^1 + T^) + 0{a') , (B2) 



32m^7r2 V 3m/ 

and higher order terms can be calculated straightforwardly. The first region that we shall 
consider when the loop momentum is small is the soft-radiation region, that we already 
found in the 1-loop diagrams contributing to the o-Ps— 87 amplitude. According to the 
hierarchy g° ~ q^/m ~ a; ~ ma we expand the propagators retaining terms up to 0{a): 

[d^q] 



-'2 



(_q2 + q^p^ _ ^) (_q2 _ g^p^ _ y) (_q2 + ^^p^ _ ^p^) 

y-ql , y-(il ^ y + '^qo^-Qo 



-q2 + qoPo - y -<f - qoPo - y -cf + go^o - ^^^o 

, 4(q-k)^ 
(-q^ + qoPo - ujPo? ' 



The integrals that result from Eq. (B3) are calculated easily with standard methods. The 
leading term in the soft-radiation region is enhanced by 1/ ^/a with respect the hard region. 
The potential region, where g° ~ q^/m y/m, also contributes when the loop momentum 
is small. The expansion of the integrand gives: 

^(P) ^ \_ [ M f^^z^ + qoPo^ 



2 P-kJ (-q2 + qoPo - y) (-q2 - qoPo - y) ^ ^ P'k 

1 Vy 



+ 0{a') . (B4) 
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The term P ■ k dominates the massive propagator depending on k, which shrinks to a 
point in the potential region. The subleading term between parentheses in Eq. (B4) arises 
from the expansion of the latter propagator, but gives a vanishing contribution because 
the numerator is proportional to one of the propagators in front. Subleading terms in 
the expansion of the propagators not depending on k generate corrections which scale as 
(go)^/(— + QqPo ~ y) ^ o?-, that we do not retain. 

The soft loop-momentum region leads to a scaleless q-integration after the q dependence 
has been expanded out from the denominators according to the scaling go q ~ "^a: 

^ I ^fA + ... = fB5l 

and the same holds for the ultrasoft region. Contributions from the soft and ultrasoft regions 
can arise in diagrams with massless propagators, like the 1-loop ladder graph of Sec. IV. 
Apart from the latter, the analysis of the threshold expansion of the 1-loop ladder diagram 
for 7^ only leads to trivial modifications to the integrations for the various regions 
found here, and can be obtained easily. 

Finally, the exact result for I2 when k"^ = can be read off the list of integrals in Ref. [6]: 

with 



L(s) = -2 (^arctany ^ j . (B7) 

One can verify that the expansion of the exact result for cu ~ ma up to terms of 0{a) is 
reproduced by the sum of the contributions of the regions in (B2),(B3) and (B4). 
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